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OBJECTIVES (CFPRC) 


Develop Verified CFD Code for Analyzing Seals 
Required Features Include: 


Applicability to a Wide Variety of Seal Confiourations 
such as: Cylindrical, Labyrinth, Face, and T?p Seals 


- Accuracy of Predicted Flow Fields and Dynamic Forces 

- Efficiency (Economy) of Numerical Solutions 

- Reliability (Verification) of Solutions 


Ease-of-Use of the Code (Documentation, Training) 
Integration with KBS 
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SCIENTIFIC CODE DEVELOPMENT 


Task 1 : Develop a 3D CFD Code (SCISEAL) for 

Cylindrical Seals 

- for Annular, Tapered, Stepped 

- Verification of Code Accuracy 

- Rotordynamlc Coefficient Calculations 


Future Tasks: Augmentation of SCISEAL 

- Incorporation of Multi-Domain Capabilities 

- Extension to Labyrinth, Damper, Face, and 
Tip Seals 


Note: 


Starting CFD Code = REFLEQS (developed by CFDRC 
under a contract from NASA MSFC/ED32) 


STATUS: 1992 WORKSHOP 

• Numerical Methods In 3D Code 

- Colocated Grids 

- High-Order Schemes 

. Rotating and Moving Grid Systems 

• Rotordynamlc Coefficient Calculation Methods 
(CFD Solutions) 

- Circular Whirl ... x 

- Moving Grid (numerical shaker) 


. Seal Specific Interface 

- Grid Generation 

- Pre-Processing 
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CURRENT STATUS 

Augmentation Effort on SCISEAL: 

• Implementation of Small Perturbation Model for flotordynamios 

' IL?! Eccentr,c a s well as Centered Seals 
Efficient, Economic Solutions 

Addition of 2-Layer Turbulence Model 


- very Small Seal Clearances -> Very Small y+ 

- 2 LauAr Mn/H Jm^ 61 lnaccurata . Low Re Model Stiffness Problems etc 
2-Layer Model Overcomes this Difficulty to Significant Extent 

Code Validation 

: LXdrfer&^.o™ nu,ar Seals ' E ~ SaaIa 

Entrance Loss Coefficients 


CURRENT CODE CAPABILITIES 


• Seals Code has: 

Finite Volume, Pressure-Based Integration Scheme 

- Colocated Variables with Strong Conservation 
Approach 

- High-Order Spatial Differencing - up to Third-Order 

- Up to Second-Order Temporal Differencing 

- Comprehensive Set of Boundary Conditions 

- Variety of Turbulence Models (k-e, Low Re k-e, 
multiple scale k-e, 2-Layer Model), Surface Roughness 
Treatment 

- Moving Grid Formulation for Arbitrary Rotor Whirl 

- Rotordynamic Coefficient Calculation Methods, CFD 
Based Centered Seals: (i) Circular Whirl (ii) Numerical 
Shaker 

■ Small Perturbation: Centered & Eccentric Seals 
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SEAL SPECIFIC CAPABILITIES 


• GUI and Preprocessor - Geared for Seals Problems 

• Easy, Quick Geometry Definition and Grid Generation 

• Four Types of Cylindrical Seals: 

- Annular, Axial Step-Down, Axial Step-Up, 
and Tapered 

• Pull-Down Menus for Problem Parameter Specification 

• One Line Commands for 

- Automatic Grid Generation 

Integrated Quantities: Rotor Loads, Torque, etc. 

- Rotordynamic Coefficients 


ROTORDYNAMIC COEFFICIENTS 


Relation Between Fluid Reaction Forces and Rotor Motion 
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Inertia (mass) 
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ROTOBDYN AMIC coefficient METHODS 

• Circular Whirl Orbit Method ' ' 


Rotor Undergoes Circular Whirl 

rm e"? £ rame Quas| - S teady Solution 
CFD Solutions at Several Whirl Frequencies 
Integration to Yield Rotor Loads 
Curve Fit to Force vs Whirl Frequency 


For Centered Rotor with Skew Symmetry Coefficient Matrices 





ROTORDYNAMIC COEFFICIENTS 
• Numerical Shaker Method 


Rotor Motion Along a Radial Direction 

Time-Dependent Solutions 

Moving Grid Algorithm for Grid Deformation 

Time-Dependent Pressure Loads -> Rotordynamlc 

Coefficients 

Can Treat Centered as well as Eccentric Seals 
Time Accurate Solutions — y Computationally Slower 
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ROTORDYNAMIC COEFFICIENT 
CALCULATIONS 


• Small Perturbation Method 

- For Centered or Eccentric/Misaligned Seals 

- Rotor Undergoes Circular Whirl with Very Small Radius 

- Resulting Perturbations In Flow Variables: 

<|> = 

- Generate Oth and 1st Order Flow Equations 

- Use Fournier Series In Time for Perturbations: 

-- Complex Form of 1 st Order Variables; 

- Flow Equations are Quasi-Steady 
. Complex Flow Perturbations Solved at Several Whirl 
Frequencies 

- Integrate Pressure Perturbations for Rotor Loads 

- Curve Fit for Rotordynamlc Coefficients 



t * nT ♦ TI2 I ■ riT ♦ 5T/B t « nT ♦ STrt 

Time-dependent solutions of the perturbation pressure 
e - 0.0, Plane at half seal length, .0.- 2.0o> 
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t » nT 


t * nT>T/8 


I - nTVTK 


t - aT*3T/B 



t-nTVra l - nT ♦5T/8 t - nT*3T/4 t = nT+77/B 

Time-dependent solutions of the perturbation pressure 
e - 07, Plane at half seal length, 0-2.5© 


2-LAYER TURBULENCE MODEL 

• Small Seal Clearances -» very Low y + Values 

• Standard Wall Functions — » Inaccurate 


• Low Re k-e Model for Very Low y + 

- can generate very stiff systems 

• 2-Layer Model Uses 

- wall functions for large y + 

- Low Re k-e model for very low y + 

• A Buffer Zone Used to Smoothly Merge the Two 
Treatments 

• Model has been Tested for a Number of Seal and 
Rotating Flow Problems 


‘Work Performed by Drs. Avva and Lai of CFDRC 


41 



SAMPLE RESULTS 

• Computation of Flow In Enclosed Rotor 
System (Dailey and Nece) 



CODE VALIDATION AND 
DEMONSTRATION 


. Code has been Validated for a Large Number of Benchmark 
Problems 

- A List of 29 Relevant Problems Included in the Interim 
Report 

• Extensive Validation Effort Conducted for Practical Seals: 

- Annular and Tapered Seals 

- Labyrinth Seals 

• Annular Incompressible Seals (Dietzen and Nordmann, 1987) 

• Long Incompressible Seals (Kanemorl & Iwatsubo, 1992) 

• Eccentric Annular Seal (Simon & Frene, 1991) 

• Annular and Tapered Gas Seal (Nelson, 1985) 

• Labyrlngth Seals Planar, (Wittig et al, 1987) 

. Labyrinth Seals, Tapered Knives; stepped 

(Tipton et al, 1986) 
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VALIDATION CASES 


1. Fully-developed flow In a pipe and channel. 


2 . 


3. 


Developing laminar flow in a narrow annulus between two cylinders 
Slug flow at inlet, fully-developed flow at outlet 

Laminar flow between rotating cylinders. Below critical Taylor number 
tangential flow only. ’ 

U^^ta‘I^ten{ nder8 ’ rOta,in0 innerCyNnder - Taylor vortex (low, 

nUmber UP *° 10 ’ 000 - C ° mpariSOnS wlth 

3-D driven cavity flow. 

Couette flow under different pressure gradients. With and without heat 
iransier. 

Planar wedge flow in a slider bearing. 

Laminar flow over a back step. Reattachment length comparison with 
experiments by Armaly and Durst. H in 

10. Laminar flow in a square duct with a 90° bend. Comparison with 
experimental data by Taylor et.al. 

1 1 . Shock reflection over a flat plate. 

12. Turbulent flow In a plane channel. Fully-developed solution at exit 
compared with experiments by Laufer. 

13. Turbulent flow induced by rotating disk in a cavity. Comparison with 

experiments by Daily and Nece. K 

14 ‘ * tator " r ° tor configuration. Comparison with 

experiments by Dibelius et.ai. 

15. Flow between stator and whirling rotor of a seal. 2-D results for 0, 0 5 and 
synchronous whirl frequencies ' ’ ana 


5. 

6 . 

7. 

8 . 
9. 
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VALIDATION CASES 


16. Flow over a bank of tubes. 

■| 7 _ Turbulent flow In an annular seal. Comparison with experiments by 
Morrison et.al. 

1 8. Turbulent flow In a 7-cavity labyrinth seal. Comparison with 
experiments by Morrison et.al. 

1 9. Turbulent compressible flow and heat transfer in turbine disk cavities 
Athavale et.al. 

20. 3-D driven cavity flow with lid clearance and axial pressure gradient. 
Control of flow through vortex Imposition. 

21. Flow in cavities on a rotor for an electrical motor. Interaction of Taylor 
vortices with driven cavity flow. 

22. Flow In infinite and finite length bearings (without cavitation). 
Comparison of calculated attitude angles with theory. 

23 Flow and rotordynamic coefficient calculation for straight, 
incompressible seals. Comparison with results from other numerical 
and analytical solutions; Dietzen and Nordmann. 

24 Flow and rotordynamic coefficients In tapered compressible flow seals. 
Comparison with bulk-flow theory results; Nelson. 

25. Rotordynamic Coefficients In a long annular Incompressible flow seal. 
Comparison with experimental data; Kanemori and Iwatsubo. 

26 Calculation of entrance loss coefficients In the entrance region of a 

generic seal. Effect of flow and geometry on the loss coefficient 

values; Athavale et.al. 

27. Flow coefficient and pressures in a 5 cavity, straight knife, look-through 
labyrinth seal. Comparison with experimental data ; Witting et.al. 

28 Flow coefficients and pressures in a 3 cavity, tapered knife, look- 
through labyrinth seal. Comparison with experimental data; T ipton 

et.al. 

29. Flow coeffients and pressures In a 2 cavity, straight-knife, stepperd 
labyyrinth seal. Comparison with experimental data; Tipton et.al. 
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LONG ANNULAR SEALS 


• Experimental Data by Kanemori & Iwatsubo (1992) 

• R = 39.656 mm, L = 240 mm, Rotor Speed = 600-3000 ram 
Clearance = 0.394 mm, Ap = 20 kPa - 900 kPa 
Specified Inlet Loss Coefficient, Ra = 1000-18000 


• Various Models Checked: 

- Whirl Method, Perturbation Method 

- Low Re k-e Model, 2-Layer Model 

- 20x15x30 grid 


DIRECT & CROSS-COUPLED STIFFNESS 



Symbols: 

Lines: 


Experimental Data by Kanemori and Iwatsubo 
Numerical Results from SCISEAL 
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DIRECT & CROSS-COUPLED DAMPING 



AP Pa 



1980 


ECCENTRIC SEAL 


• Experimental Data by Falco et al (1984) 

Numerical Data by Nordmann (1987), Simon & Frene (1991) 


• Radius = 80 mm, Length = 40 mm, e = 0. 1-> 0.7 
4000 rpm, Ap = IMPa, Entrance Loss Coefficient = 0.5 


• Physical Models 

- Standard k-e Model 

- Small Perturbation Method 

■ 12x6x30 grid 


DIRECT STIFFNESS 
COEFFICIENT, K w , K 77 
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CROSS-COUPLED STIFFNESS, 

: YZ' 


K\/7, K^y 


3. OEfOSn 



QO 


K 


yz 



5. Ef06" 
N/m 
4.Et06 

3. Ef06 

ZEf06 


Nelson and Nguyen __J 

Simon and Frene 1 * 

Nordmann 


— 

02 


— t — 

04 


Q6 

Eccentricity e 


Q8 


oo 


K* 


QO 


Q2 


~r 

Of 


04 

Eccentricity e 


DIRECT DAMPING, C vv , C a 




Eccentricity, e Eccentricity, e 



49 


STRAIGHT LABYRINTH SEAL 


. Experimental Data by Wittig et al (1987) 

• 5 Cavity, Planar Look Through Seal 

• Physical Models 

- 30x30 Cells In each Cavity, 8/1 2 Cells In Gap 

- Compressible Flow, Standard k-e Model 

- Specified Pressure Ratio Across Seal 

• Results: Numerical Results Compared with 
Experimental Data 

- Pressure Along the Seal Length for Different Tip 

- Mass Flow Rates at Different Tip Gaps and Pressure 
Ratios 


RESULTS FOR STRAIGHT LABYRINTH 
SEAL . 


Pressure 


Mass Flow 



i 
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Pstator^inlet 


STEPPED LABYRINTH SEAL 

• Experimental Data by Tipton et al (1 986) 

• 2 Cavity, Planar, Stepped Labyrinth Seal 

• Physical Models: 

Compressible Flow, Specified Pressure Ratio 

- Standard k-e Model 

- 26x53, 26x62 Cells In Cavities; 1 0 Cells In Tip Gap 

• Results: Numerical Results Compares with 
Experimental Data 

- Pressure Along Stator and Rotor Surfaces at One 
Pressure Ratio 

- Mass Flow Rates at Different Pressure Ratios 
‘Work Performed by Dr. Makhijanl of CFDRC 


RESULTS FOR STEPPED LABYRINTH 
SEAL 
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• Measure of Flow Losses at Entrance Region 
. SCISEAL Used to Compute £, with CFD Solution 


. Variation of £ with 

- Axial Reynolds Number 

Seal Clearance-to-Radlus Ratio 

- Entrance Gap-to-Clearance Ratio 




Flow 


entrance 

region 


I 

seal 


. Physical Models 

- Incompressible Flow, Standard k-e Model 

- Fully Developed Flow Upstream, Pressure 

50 Cells In Axial Direction, 5 In Clearance, 30 or 50 In 
Entrance Region 


FI OW GEOMETRY FOR ENTRANCE LOSS 
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RESULTS 


Table 1. Entrance Loss Coefficients, Radius / Clearance - 50 


Entran 

ce Gap/Clearance * 50 

I Entrance Gap/Clearance - 1 00 1 

Ugx nVs 

Re ax 

$ 

Uax nrVs 

Re ax 

$ 

10.814 

16.232 

21.619 

26.942 

10377 

15484 

20746 

25854 

0.471 

0.431 

0.414 

0.406 

10.82 

10.24 

21.66 

27.06 

10384 

15584 

20785 

25970 

0.490 

0.488 

0.482 

0.48 

Table 2. Entrance Loss Coefficients, Radius / Clearance = 

PA-ti2-34 Q 
= 100 

Entranc 

)e Gap'Clearance - 50 

Entrance Gap/Clearance * 1 00 1 

Uax m/s 

Re ax 

5 

Uax 

Reax 

5 

10.80 

16.56 

21.595 

26.67 

32.27 

43.062 

5181 

7945 

10361 

12796 

15484 

20667 

0.562 

0.54 

0.526 

0.51 

0.493 

0.478 

10.797 

16.176 

21.55 

26.934 

32.24 

42.533 

5167 

7761 

10339 

12664 

15469 

20408 

i 


PA-aa^o 


Table 3. Entrance Loss Coefficients, Radius /Clearance = 150 


Entrarx 

:e Gap/Clearance « 50 

I Entrance Gao/Clearance - 100 1 

Uax 01/8 

^ e ax 

* 

Uax m/s 

Re ax 

? 

10.82 

3461 

0.66 

10.75 

3438 

0.68 

16.19 

5178 

0.65 

16.09 

5146 

0.66 

21.49 

6874 

0.647 

21.47 

6874 

0.65 

26.74 

8553 

0.637 

26.81 

8553 

0.648 

32.25 

10315 

0.628 

32.176 

1 10292 

0.64 

48.33 

15461 


47.87 

15315 

0.63 

64.487 

20630 

0.595 | 

64.165 

20630 

0.624 


RELATED CFD RESULTS 


REFLEQS-3D Used for Rotating Flows 

Pump C^ns^Xn?) Cen,r,,u 9 al lm P eller < For “SFC 
Techniques 0 & SC,SEAL have Slml,ar Numerical 


• SCISEAL In Narrow, Long Channels 


Suitable for Cooling Channels In Rocket Nozzles 
Heat Transfer & Flow Calculations 
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Secondary flow patterns In a long, narrow channel 

(Velocity vector *ae and croat-aaction size* not to acala) 


CONCLUDING REMARKS 

• A 3D CFD Code, SCISEAL, Being Developed and Validated 

Current Capabilities Include Cylindrical Seals 

• State-of-the-Art Numerical Methods 

- Colocated Grids 
High-Order Differencing 

- Turbulence Models, Wall Roughness (In progress) 

• Seal Specific Capabilities 

Rotor Loads, Torques, etc 

• Rotordynamlce Coefficient Calculations 

Full CFD Based Solutions - Centered Seals 
Small Perturbations Method - Eccentric Seals 

• Extensive Validation Effort 
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WORK PLANS FOR NEXT YEAR 

• Consolidate Current Models 

• Include Multi-Domain Solution Methodology 

- Efficient Solutions for Complicated Flow Geometries 

- entrance region & seal clearance 

- stepped and straight labyrinth seals 

- face seals 
-- tip seals 

- conjugate heat transfer 

- Increases Code Flexibility 

- Technology Already Developed but Requires j 

Adaptation and Testing for Seals 

4 

• Continue Work on Labyrinth Seals 

• Validation/Demonstration for Practical Seal j 

Configurations j 



Single Domain Grid 



Multidomain Grid 
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Stepped Labyrinth Seal grids 



♦ 


Single Domain grid 



Multidomain pid 
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